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Abstract 

We study the linearized hydrodynamics of a two-component fluid mem- 
brane near a repulsive wall, via a model which incorporates curvature- con- 
centration coupling as well as hydrodynamic interactions. This model is a 
simplified version of a recently proposed one [J.-B. Manneville et al. Phys. 
Rev. E 64, 021908 (2001)] for non-equilibrium force-centres embedded in fluid 
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membranes, such as light-activated bacteriorhodopsin pumps incorporated in 
phospholipid egg phosphatidyl choline (EPC) bilayers. The pump/membrane 
system is modeled as an impermeable, two-component bilayer fluid membrane 
in the presence of an ambient solvent, in which one component, representing 
active pumps, is described in terms of force dipoles displaced with respect 
to the bilayer midpoint. We first discuss the case in which such pumps are 
rendered inactive, computing the mode structure in the bulk as well as the 
modification of hydrodynamic properties by the presence of a nearby wall. 
These results should apply, more generally, to equilibrium fluid membranes 
comprised of two components, in which the effects of curvature-concentration 
coupling are significant, above the threshold for phase separation. We then 
discuss the fluctuations and mode structure in steady state of active two- 
component membranes near a repulsive wall. We find that proximity to the 
wall smoothens membrane height fluctuations in the stable regime, resulting in 
a logarithmic scaling of the roughness even for initially tensionless membranes. 
This explicitly non-equilibrium result is a consequence of the incorporation of 
curvature-concentration coupling in our hydrodynamic treatment. This result 
also indicates that earlier scaling arguments which obtained an increase in the 
roughness of active membranes near repulsive walls upon neglecting the role 
played by such couplings, may need to be reevaluated. 
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I. INTRODUCTION 



Amphiphilic molecules, in polar solvents such as water, form symmetric bilayer membrane 
phases at sufficient concentration [1-4]. Such membranes are typically fluid in nature, with 
short-ranged positional order, although more complex forms of ordering are possible [5-7]. 
The equilibrium conformations of such single-component membranes are governed by the 
energy cost for bending the bilayer, provided the surface tension can be neglected, as is the 
case if the membrane is self-assembled [4,8,9]. Our understanding of the static and dynamic 
properties of such single-component membrane systems in equilibrium is now fairly detailed 
[3,10-17]. 

More complex multi-component membranes have also been synthesized and studied 
[9,18-20]. Such membranes are especially important in the biological context, since cell 
membranes are often usefully idealized as a bilayer fluid-mosaic comprised of over a hun- 
dred different types of lipid and protein constituents [21,22]. Biological membranes, however, 
often have additional complicating features: An underlying network of cross-linked proteins 
anchored to the bilayer and associated with the cell cytoskeleton typically lends such mem- 
branes a non- vanishing (although small) shear modulus [22]. Another aspect, specific to the 
biological context, is the presence of non-equilibrium driving forces: An important class of 
trans- membrane proteins are the "ion pumps" , molecules which consume energy derived from 
ATP hydrolysis, electro-chemical gradients or light, and undergo conformational changes. 
Such pumps maintain electro-osmotic potential gradients across the cell membrane by con- 
trolling the flow of ions such as K + and Na + . Since energy must be supplied externally, the 
driving of the pump is a process which occurs out of thermal equilibrium. Thus, the situa- 
tion of active pumps diffusing in a fluid membrane matrix is an intrinsically non-equilibrium 
problem whose behavior represents a novel class of non-equilibrium steady states. 

The appropriate statistical description of active or "driven" biological membrane systems 
has attracted recent experimental and theoretical attention [23]. It has been realized that 
non-equilibrium behavior may underlie aspects of bio- membrane dynamics [24,25] previously 
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attributed purely to equilibrium thermal fluctuations, such as the "flicker" phenomenon in 
erythrocytes [26]. Recent micropipette experiments on light-activated bacteriorhodopsin 
(BR) pumps incorporated in phospholipid (EPC) bilayers find that area fluctuations in 
such active membrane systems can be phenomenologically described in terms of an effective 
temperature which exceeds the true physical temperature by up to a factor of 2 [27,28]. 
In these systems, the transition between active and inactive states can be easily (and re- 
versibly) tuned, providing a remarkable window into the dynamics of steady states away 
from equilibrium [28-31]. 

Analytical calculations which deal with the complex nature of a typical biological mem- 
brane in toto are difficult, if not impossible. We may, however, hope to gain useful intuition 
by working with simpler models. We idealize a typical biological membrane here as com- 
posed of principally two different types of molecular constituents, the lipids which constitute 
the bulk of the bilayer membrane and proteins which diffuse freely on the membrane surface. 
Such simple lipid/pump protein systems can be reconstituted and studied in vitro, as in the 
micropipette experiments referred to above. Our results should apply, most directly, to such 
experiments. 

The question we address in this paper is the following: How are the fluctuations of 
such a two-component membrane [32], both in equilibrium and out of equilibrium, affected 
by the presence of a nearby repulsive wall? Fluctuations of a single-component membrane 
bounded on one side by a repulsive wall have been considered in earlier work [3,33]. More 
recently, a similar problem for active membranes has been studied in an influential paper by 
Prost, Manneville and Bruinsma [29,30], who suggest that steady-state fluctuations of non- 
equilibrium membranes near repulsive walls can be amplified substantially relative to the 
equilibrium case. This amplification is studied via a scaling treatment which incorporates 
non-linear effects. Our starting model differs from theirs in several respects, in particular 
in the way we describe activity, in our incorporation of curvature-concentration coupling, 
as well as in our evaluation of the relative importance of permeative and hydrodynamic 
effects at the length scales accessed by typical experimental probes. We study the linearized 
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hydrodynamics of such active pump-membrane systems in this paper, ignoring the role of 
non-linear effects. We will describe the similarities as well as the differences in the results 
we obtain. 

We answer the question posed in the previous paragraph first through a calculation of the 
correlations and linearized mode structure in thermal equilibrium of a model impermeable 
two-component fluid membrane, incorporating hydrodynamic interactions. This amounts, 
in terms of the model for biological membranes discussed above, to assuming that the pump 
proteins are "passive" constituents of the membrane. While motivated here in a biolog- 
ical context, such models are also useful outside this specific context. Physics similar to 
that described in the case where the pumps are rendered passive should be generic to two- 
component systems composed of (i) a "symmetric" (lipid) constituent forming an up-down 
symmetric bilayer phase in isolation, and (ii) an asymmetric (protein) constituent, in which 
either the shape of the molecule or its location with respect to the bilayer favors a local 
spontaneous curvature for the bilayer. For concreteness, we will often use the terms lipid 
and protein when referring to constituents of types (i) and (ii), following our earlier dis- 
cussion, but stress that the results obtained here in the equilibrium case are also applicable 
more generally. 

For protein molecules which lack up-down symmetry, we may distinguish between 
molecules oriented parallel to the local normal and those which are oriented antiparallel 
to the normal. We label these proteins as "+" and "-" for convenience; the relatively small 
flip-flop rate for protein transfer across leaves of the bilayer ensures that these labels are con- 
served across short and intermediate time scales. The difference between the local densities 
of + and - proteins defines a "signed " protein density field tp. To describe the incorpo- 
ration of such asymmetric proteins into the lipid bilayer, we model the protein as a rigid 
rod of length w 1 ' + wK The protein is taken to situate itself asymmetrically with respect 
to the bilayer midpoint; a section of the protein lies on one side of the bilayer midpoint 
while a section lies on the other side. This relatively simple model of the asymmetry is 
convenient for analytic computation. 
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To lowest order, the two-component character of the membrane manifests itself via the 
existence of a curvature-concentration coupling - fluctuations in ip influence the mean curva- 
ture in that region (see Fig. 1). We will work with balanced membranes for which (ip) = 0. 
For inactive pumps, corresponding to the equilibrium case, the curvature-concentration cou- 
pling term in the free energy effectively shifts the bending rigidity of the membrane from k 
to n e = k — (kH) /xo 1 - Here H is the coefficient coupling curvature to concentration and 
Xo 1 an inverse compressibility [31]. 

In closed vesicular structures found in biological contexts, while the bilayer constituting 
the vesicle may be symmetric, chemical and physical environments within and outside the 
vesicle can differ considerably. Thus, a pump protein can, in principle, distinguish the side of 
the bilayer exposed to the inner volume from the side which is exposed to the outer volume. 
For a large enough vesicle, flat on the large length scales of relevance to a hydrodynamic 
calculation, the pump molecule can thus tell up from down, where up and down are defined 
arbitrarily with respect to the average position and orientation of the bilayer midpoint. As 
an extreme limit, we may consider the case in which the pump molecules preferentially insert 
always on the same (say up) side, so that a section of the protein molecule always lies 
above the bilayer midpoint while a section lies below. 

A recent study of active membranes works with this model of the architecture of a 
membrane-pump system, modeling activity in terms of dipole forces associated with the 
ends of the pump, as we detail below [28]. The associated forces may point either outward 
at the ends (up pumps) or inward at the ends (down pumps); see Fig. 2. In this model, 
up and down pumps are physically distinguishable. The main results presented here are 
appropriate to this model. 

Note that this distinction between "above" and "below" is artificial for strictly symmetric 
bilayer membranes, in which pump proteins may insert on either leaf of the bilayer with 
equal probability. In an Appendix, we discuss results for an alternative model, in which the 
protein attachment respects the "up-down" symmetry of the bilayer; the center of mass of 
the protein can be displaced either above or below the bilayer mid-plane as shown in Fig. 3. 
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The head and tail distances from the bilayer midpoint are fixed at wi and w s respectively. In 
this model, up pumps and down pumps are physically indistinguishable; their nomenclature 
follows from the side of the bilayer on which they are inserted. 

We will present results for non- equilibrium membranes in the bulk, in which case the 
pumps are active, and for the case in which fluctuations of such an active membrane are 
bounded on one side by a repulsive wall. This study complements results by Prost, Man- 
neville and Bruinsma for a model in which curvature-concentration coupling is ignored, an 
effect explicitly incorporated in later work [31]. In the active case, the pumps exert forces 
on the membrane and the surrounding fluid. Since these forces are internal forces, they 
must cancel when integrated over the size of the pump [34]. Following Manneville et al. in 
Ref. [28], we consider the simplest model for pumps consistent with this requirement: the 
pumps are taken to be dipole force centres with positive and negative force centres located 
asymmetrically with respect to the midpoint of the bilayer. 

Our other major assumptions are the following: We assume that our observations are 
conducted on scales such that the membrane may be considered to be impermeable to the 
solvent in which it is embedded. This assumption follows from the work of Manneville 
et al. who observe that for active terms permeative effects can be neglected vis. a vis 
hydrodynamic ones over a substantial range of length scales upto microscopic ones. We 
derive our results under the assumption that the active membrane is impermeable and 
extend this assumption to the case when the membrane is in the passive state, to smoothly 
interpolate between results in the two regimes. (In the Rouse or free draining limit, where the 
membrane dynamics decouples from that of the fluid, the membrane is implicitly permeable). 
We take the solvent to be incompressible and work at small Reynolds number, neglecting 
the inertia of the fluid. We use linearized equations of motion for the hydrodynamic velocity 
field, corresponding to the Stokes limit of the Navier-Stokes equations. 

We now summarize our principal results. For equilibrium, impermeable, two-component 
membranes bounded on one side by a wall placed at a distance d from the membrane plane, 
we obtain the following modified mode structure to leading order, 
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(1.1) 



In our notation, D is the diffusion coefficient of the proteins on the membrane, a the surface 
tension of the membrane, k the bending rigidity of the membrane and rj the viscosity of the 
fluid. Here, 



is a "crossover function", which describes how membrane fluctuations interpolate between a 
regime where physical interactions between the membrane and the wall cannot be neglected, 
and a regime where such fluctuations do not feel the presence of the wall [33]. In the 
A;_|_ ^> 1/d limit, the above result reduces to one derivable in the bulk case, as is physically 
sensible. 

To study non- equilibrium membranes we model active pumps, following Manneville et al., 
as "asymmetric" force dipoles with force-centres located at distances above the bilayer 
midpoint and w^- below the bilayer midpoint. These distances are, in general, unequal. The 
magnitude of the force exerted at each of these force-centres is /. The pumps diffuse on the 
membrane and are characterized by a mobility A. 

In the presence of a wall, for non-equilibrium membranes, it has been argued that static 
height fluctuations are amplified and an interesting speculation is that such an amplification 
may, in fact, speed up the process of membrane fusion when two non- equilibrium membranes 
are brought close together [29]. When the wall is at a distance d from the membrane, our 
calculation yields (for a tensionless membrane) 






Xo 



(1.3) 



where Q aw is the crossover function 




(1.4) 
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fl(k±,w\w^) is proportional to the degree of asymmetry (w^ — w^). In a regime where 
fluctuations feel the wall, non-equilibrium effects do not depend on the asymmetry factor 
but only on u^, as seen above. Here again, curvature-concentration coupling plays a crucial 
role in the manifestation of non-equilibrium effects in the linearized mode structure. 

We find that in the presence of the wall, the first term in oj\ is the most important since 
it is of order 0{k\). The existence of this mode depends crucially on the presence of both 
activity and curvature-concentration coupling. It will be shown later in this paper that this 
mode causes active membranes near walls to relax faster than their equilibrium counterparts 
and causes a smoothening of membrane fluctuations; in the presence of a wall, the roughness 
of tensionless, active membranes averaged over a patch of size LxL, ((/i 2 (L))), scales as ln(L) 
instead of L 2 . Thus the wall acts to smoothen active membrane fluctuations in a tensionless 
membrane by inducing an apparent surface tension. This result is a purely non-equilibrium 
one, since the apparent surface tension vanishes when the strength of active forces is set to 
zero. It also relies crucially on the existence of a curvature-concentration coupling, since 
the apparent surface tension also vanishes when the curvature and concentration fields are 
decoupled. We show that this effect can be understood in terms of the dynamics of the fast 
relaxing tp field which drives the height field h(r±, t) via activity and curvature-concentration 
coupling. 

The outline of this paper is as follows: Section II describes our basic coarse-grained free 
energy and reviews the calculation of the static properties of a two-component membrane 
with curvature-concentration coupling. In Section III (a), the free-draining approximation 
(Rouse model) is used to discuss the dynamics. In Section III (b), we incorporate the 
effects of the hydrodynamic velocity field of the solvent. Section IV describes the mode 
structure in the presence of a wall for the equilibrium problem. In Section V, we discuss the 
mode structure and correlation functions of active two-component membranes. Section VI 
discusses how these properties are modified by the presence of a wall. In Section VII, we 
summarize the conclusions of this paper. An Appendix discusses the details of an alternative 
model for active two-component membranes, with up-down symmetry. 
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II. STATICS 



We model the membrane as an infinitesimally thin two-dimensional surface embedded in 
three-dimensional space. Points on this surface are specified by a three-dimensional vector 
R(w), with components Ri(u), i — 1 . . . 3. We use Monge gauge, valid for nearly flat surfaces, 
in which the surface is specified by its height above a two-dimensional flat plane. This plane 
is parameterized via Euclidean coordinates i.e. u = (x,y) = r_i_; R(rj_) = (x,y,h(x,y)). In 
this gauge, the curvature is H = —V 2 h, to lowest order. 

The number densities of the lipid molecules, + proteins and - proteins are specified 
through fields pg, p + {u,i) and p~(u,t) respectively. These are the physical densities of the 
lipids and proteins on the membrane. The projections of these densities on the x — y plane 
are obtained by multiplying the physical density fields by the metric factor \Jl + (Vh) 2 . 
To lowest order in gradients of h, physical densities and projected densities coincide. The 
scalar fields (f>(u,t) and ip(u,t) are defined by <j)(u,t) = p + {u,i) + p~(u,t) and i/)(u,t) = 
p + (u,t) -p~(u,t). 

The free energy for the membrane-protein system consists of a part arising from protein 
and lipid densities (the matter part) and another which arises as a consequence of membrane 
fluctuations about the flat state. Given equilibrium concentrations of lipid and protein (p ), 
(<f>(u,t)) = (p + (u,t)+p~(u,t)) = 0o , {ij>(u,t)) = (p + (u,t)-p-(u,t)) = ip and fixing p ~ p 
and ~ 0o, the appropriate gauge- invariant free energy of the membrane with asymmetric 
proteins is 

F = J d 2 u^[a + l -x*\i> ~ ^o) 2 ] + \kJ d 2 u^(H - H ) 2 . (2.1) 

where a is the surface tension and k is the bending rigidity. To lowest order we take 
H = Hip. To leading order in gradients of h, we have 

F=\j d 2 r ± [ K (V 2 h(r ± , t)f + a(Wh(r ± , t)) 2 - 2Hki/>(t ± , t)V 2 h(r ± , t) + X q V(r±, t)\ 

(2.2) 

10 



yielding a local force density of f(rj_) = —5F/5h . The local force density is evaluated at 

V 

constant density [36]. 

The partition function is [41] 

Z = J 'v[h]V[iP]e- pF ^\ (2.3) 

which yields, on integrating out the ip field, an effective free energy for the height field 

FeffH 

Feff[h] = \ + Kk{)\h(k ± )\ 2 - { -^kl\h(k ± )\ 2 . (2.4) 

2 k ± *o 

The equipartition theorem yields the equilibrium correlation function for the height field 
fluctuations (about the equilibrium flat, "mixed" phase) as 

<"*( k J^» - K^gj ' (2 ' 5) 

where K e = k—(kH) /xo _1 - The asymmetric components act to reduce the rigidity modulus 
from k to K e as a consequence of curvature-concentration coupling [31], an effect easily 
understood on physical grounds. We assume that H is small enough to keep K e positive, in 
which case the membrane is linearly stable. Instabilities arising at equilibrium due to large 
values of H have been discussed in [40] . 

III. MODE STRUCTURE IN THE BULK 

We now study the dynamics of a two-component membrane with curvature-concentration 
coupling and suspended in a solvent. We will first present results in the Rouse model limit 
in which the background fluid drains freely through the membrane and effects due to the 
hydrodynamics of the background fluid velocity field can be ignored [43]. 

A. Rouse Model 



In this model the fluid velocity field decouples from the membrane height field and ip. 
The membrane is implicitly permeable and the height field fluctuations are governed by a 
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balance between the elastic forces and local (permeative) friction. An equation of the form 



SF 

d t h(r ± ,t) = -r— + f m , 



(3.1) 



follows, where T is a kinetic coefficient and f m is a thermal equilibrium noise whose statistics 
ensures that the static height-height correlation function is Eq. (2.5). 

The proteins diffuse freely on the surface of the membrane, Since ip is locally conserved, 
it obeys an equation of continuity of the form 



SF 

<9^(r ± ,t) = AV 2 — + V-f^, 



(3.2) 



with A = D/xo^ 1 , D is the diffusion coefficient of the proteins on the membrane. A contri- 
bution due to the in-plane current of the pumps has been dropped; this is typically small and 
can be neglected to linear order. This term contributes if we go beyond the strictly linear 
treatment by replacing ip 2 by (ip 2 ), yielding a convective term which gives rise to waves whose 
speed is independent of the wavevector [23,31]. The last term is a conserving Gaussian noise 
with (f^,(rj_,t)) = and correlations (f^i(r±,t)f^j(r' ± ,t')) = 2Ak B T5ij5(rj_ — r'j_)5(t — t'). 
Equations (3.1) and (3.2) can be written as 



^ h(k±,uj) ^ 



AKHk 4 ± Axo lk l 



\ 



J 



( 



h(k ± ,u) 



\ ( 



+ 



fn 



\ 



(3.3) 



Solving these coupled equations yields 
(h*(k ± ,cu)h(k ± ,u;)) = ^ 



(fmfm)(uJ 2 + T7 12 ) + 2(TKHfAk B Tkl 



{-^ + r fc " V - AT(KHfklY + u *(t? + O 2 ' 

where r^ 1 = T(ak 2 ± + nk\) and r^ 1 = A\o 1 k\. 

On integrating over u>, this correlation function should reduce to the same form as Eq. 
(2.5). This constrains the equilibrium thermal noise correlations to be the following: 



(3.4) 



(fm(r ± ,t)f m (r' ± ,t')) = 2Tk B T5(r ± - r' ± )5(t - t'). 



(3.5) 



In addition, we impose (/ m (r_|_,t)) = 0. 
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To obtain the mode structure, we drop the noise terms in Eq. (3.3) and find the eigen- 
values of the matrix. This yields 

u u U2 = ^ (r(<7ifel + Kki) + Dkl) ± \ (V^V^I + Kki) + Dk lf - ±DTk\{ok\ + /c e fci)) 

In the absence of curvature-concentration coupling these modes decouple and are 

= T(ak\ + Kki), Dhj_, the dispersive modes of a tense membrane [35]. For tensionless 

membranes, the presence of curvature-concentration coupling modifies the mode structure 

to 

U! « TK e ki, uj 2 « DA; 2 + ^4r rA; l- ( 3 - 6 ) 



B. Zimm Model 

Next, we examine the effect of solvent hydrodynamics on the dynamical correlation func- 
tions of the protein density field and the membrane height field [43]. Since no external forces 
act on the (fluid+membrane+proteins) system, momentum conservation requires [42,44] 

d t 9i = -VjiTij, (3.7) 

where g is the momentum density and -k^ is the stress tensor. The stress tensor takes the 
form TTij = P5ij + pViVj — r](\7 iVj + V jVi) + Sij, where v(r,t) is the fluid velocity, P(r,t) is the 
pressure field and «Sy describes the contribution to the stress tensor arising from membrane 
conformations (the V 2 in the above equation is evaluated in 3 dimensions). In the Stokes 
limit, the non-linear term can be dropped. At low Reynolds number, we neglect the inertial 
term on the left hand side of Eq. (3.7) leading to V 'j-Kij = 0, and thus the equation governing 
the hydrodynamic velocity field 

r]V\(r, t) - VP(r, t) - F(r, t) + f h (r, t) = 0. (3.8) 

F(r,t) is the force acting on the fluid due to membrane stresses Sy and f h (r, t) is the 
equilibrium thermal noise present in the bulk of the fluid. This noise obeys (fh(r, t)) = 
[28,29] and is related to the viscosity via the fluctuation-dissipation theorem. 
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(f hi (r,t)f hj (r',t')) = 2k B Tr 1 {-5 lJ V 2 + 8^)5^ - r')S{t - t'), (3.9) 

where (i, j) index the components of the force. Eqns. (3.8) are to be solved together with 
the incompressibility condition 

V-v(r,i) = 0, (3.10) 

to obtain the mode structure of an impermeable two-component membrane in the presence 
of an incompressible solvent. 

We employ a useful convention due to Seifert in which a basis spanned by the z- 
component and the longitudinal and transverse in-plane components of the velocity is used 
[3]. (This convention will be useful further on in this paper when we discuss membrane 
fluctuations in the presence of a confining wall). We solve the above equations to obtain the 
transverse (t), longitudinal (I) and z components of the fluid velocity in terms of the t,l and 
z components of the force F acting on the fluid [38] . The t and I components of the velocity 
lie in the plane of the membrane and are perpendicular and parallel to k^, the wavevector 
describing fluctuations of the membrane. Some algebra yields these components as [3] 

v z (k ± ,z) = -±- [°° dz'e' k ^'- z \[(l + k ± \z-z'\)F z (k ± ,z') + tk ± (z'-z)F l (k ± ,z% 

Vi(k ± ,z) = r dz'e- k ^'- z \[{l - k ± \z - z'DF^z^+tk^z' - z)F z (k ±: z')} (3.11) 

Ar/k^ J-oo 

T*(kj_, z) = -—\ dz'e- k ^'^2F t (k ± , z'). (3.12) 

Note that the transverse components can be ignored since they do not couple to the other 
components. We can now calculate f 2 (kj_,0) by supplementing the above equations with 
boundary conditions. 

In the absence of permeation the membrane is advected by the fluid. The velocity of any 
point on the membrane is the ^-component of the velocity of the fluid at that point [28]. 
The equation of motion for the height field of an impermeable membrane is therefore 

dth(T ± ,t)=v g (r ± ,0,t), (3.13) 
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where v z (r±,0,t) is the velocity of the surrounding fluid at the mean position of the mem- 
brane. The force acting on the fluid is of the form 



F(k ± ,z)=5(z)[-"—z + Frk ±i 



6F , 
'Jh' 



(3.14) 



where F[" is the longitudinal component of the force acting on the fluid due to the membrane. 
The elastic force density due to the membrane is —SF/5hz . We solve for v z (k±,0) using 
Eqns. (3.11), together with the boundary condition that the membrane is incompressible 
[36] i.e. the in-plane divergence of the velocity field vanishes at the membrane. This implies 
that vi(k.±, 0) = 0. This boundary condition gives F™ = 0, yielding 

1 5F 



v x Qc±,0) 



(3.15) 



Ar/k± 5h 

The bulk thermal noise present in the fluid also contributes to v z (k±,0). This contribution 
is denoted by f z , with 



(3.16) 



2nr]k 2K n 

We can use the noise correlation function given in Eq. (3.9) to show that the noise f z has 
zero mean and that its variance in Fourier space is 2k B T '/ '4r)k±. 
Finally, Eq. (3.13) becomes 

1 5F 



d t h 



+ fz(k±,o,t). 



(3.17) 



Ar]k ± 5h 

The density difference field obeys Eq. (3.2). Eq. (3.13) and Eq. (3.2) can then be written as 



1 h(k ± ,u;)\ 



I 



V 



Solving these equations gives 



\ 








1 fzfa,^ 








+ 




/ 




M ) 




ik±.i^ ^ 



(3.18) 



{2k B T/4r]k ± )(uj 2 + tI 12 ) + (k B TA(KHy/8rj 2 )k 



2U4 



(h*(k ± ,u)h(k ± ,u)) = 



(-" 2 + r h \ l - (A(KHf/A V )kl) + ^( T -i + r ^) 2 
The parameters appearing in the above equations are 



(3.19) 
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! -(akl + Kk\),T7 l =K x - l kl. (3.20) 



Integrating this result over uj recovers for us the equal time correlation function for height 
fluctuations about the equilibrium, flat, "mixed" phase 

(h*(k ± , t)h(k ± , t)) = {ak 2*f Kek Ay (3-21) 

as expected. The roughness of the height field at macroscopic length scales L, defined 
through (h 2 (L)) = fyl d 2 k±(h(k±)h*(k±)) scales as 

f hBT L 2 if a = o, 
(h (L)) = { Ke (3.22) 
( ^HL/a) if a ? 0, 

where a is a microscopic cut-off length. 
The mode structure is then 

U U U2 = — (<rk± + kAm. + 4-qDk 2 ,) ± — ( J (ak ± + kB + lr)Dk\f - 16r]DUj£ + ak 3 ± )) ■ 
S// v 8// \ ' ) 

In the absence of curvature-concentration coupling, = (ak± + Kk\)/Arj, Dkj_, i.e. the 

decoupled modes of the membrane height field and the protein concentration field in the 
presence of hydrodynamic correlations. For tensionless membranes with non-zero curvature- 
concentration coupling, the modes are 

Ul ~ ^ U2 ~ D kl + ^S-kl. (3.23) 
At] 4^xo _1 



IV. MODE STRUCTURE IN THE PRESENCE OF A WALL 

Consider a fluid bilayer membrane with shape asymmetric constituents, suspended in a 
fluid bounded on one side by a wall. Our model for the membrane-pump system follows that 
of Manneville et ai. In our co-ordinates, the membrane fluctuates about the z = d plane 
and the wall is at z = 0. The membrane is kept at distance d from the wall by a repulsive 
interaction with the wall and a suitable external pressure which prevents it from floating 
away into the bulk. 

16 



The relaxation rate of an equilibrium membrane is reduced near a wall; in tandem, the 
equilibrium noise fluctuations are required to have a reduced amplitude as a consequence 
of the fluctuation-dissipation theorem [29]. On short-length scales, however, membrane 
fluctuations do not see the wall and the height-height correlation function is identical to 
that in the bulk. This can be used to determine the form of thermal correlation functions 
in the presence of the wall. 

We work in the linear regime where (h 2 (L c )) d 2 , defining a collision length L c . Beyond 
this length the membrane senses the presence of the wall and the non-linearities neglected 
in our treatment become important; 1/L C is the smallest wave vector above which the 
linearization condition \h(r±,t)\ <C d is valid. From Eq. (3.22) we can estimate the collision 
length of an equilibrium membrane near a wall to be 



L c (d) 



x P^d if a = 0, 

a exp(ad 2 /k B T) if a ^ 0. 



We now solve Eq. (3.8) to obtain the velocity of the fluid at the position of the membrane. 
The membrane and wall are introduced in Eq. (3.8) as external forces acting on the fluid, 
imposing the required boundary conditions. The incompressibility of the fluid is used to 
solve for the components of the velocity in terms of the components of the external force. 

We express the components of the velocity in the "mixed" basis of Seifert [3] as before. 
The force exerted by the membrane and wall on the fluid is of the form [3] 

F(k ± , z) = 5(z - d)[- 5 -^z + F/"kJ + 5(z)[F?z + ifkj, (4.2) 

where F™ and F™ are the z and longitudinal components of the force exerted by the wall 
on the fluid at z — 0. The quantities F™, F™ and F™ can be evaluated using the boundary 
conditions: all components of the velocity should vanish at the wall (no-slip boundary 
condition) [44]. This imposes 

v z (k ± , 0, t) = u,(k_L, 0, t) = v t (k ± , 0, t) = 0. (4.3) 

These conditions, together with membrane incompressibility (vi(\a±,d,t) = 0), provide 
the boundary conditions for the velocity field of the fluid at the membrane. This gives 
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v z (k±, d, t) = —Q%- The equation of motion for the height field is then obtained from Eqn. 



(3.13) 



where 



G(k±, kj_d) 



5F ~ 

d t h(k ± ,t) = -g— + f z (k ± ,d,t), 



2(sinh 2 (A; ± d) - (k ± df 



(4.4) 



(4.5) 



4:T]k± \smh. 2 (k±d) - (k±d) 2 + kj_d + sinh(k±d) cosh(k±d) J ' 
d 3 k 2 J12 V if k± < l/d t, 
l/Ar]k ± if k± > l/d. 

We note that the term marked f is incorrect for arbitrarily small wave vectors. The relevant 
cut-off is the inverse of the collision length L c since linearized calculations are not valid at 
wave-vectors smaller than this value. 

The contribution to the velocity field at the membrane due to the noise in the fluid in 
the presence of the wall is f z . Following the procedure used in the derivation of Eq. (3.18), 
we get 



1 h(k ± ,u)\ 



4>(k ± ,u) 



( Q{ak\ + Kk A ± ) G{KHkl) ^ ( 



V 



AnHki AXo^l j 



h(k ± ,u) 
ip(k ± ,uj) 



fz(k ± ,u] 



ik i .f, 



(4.6) 



V> j 



This matrix equation is solved to yield 



(h*(k ± ,u;)h(k ± ,u;)) 



(4.7) 



where t^ 1w = Q{ak\ + K,h±). 

The hydrodynamics of the solvent cannot change equilibrium correlation functions. Thus, 
demanding that equal time correlations of the h field are the same as in Eq. (3.21), we obtain 
the noise correlation function of f z as 



(/,(k ± ,f)) = 0, 
(f z (k ± ,t)f z (k' ± ,t')) = 2k B TG(k ± , k ± d)5(k ± + k' ± )S(t - t'). 



(4.8) 
(4.9) 
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Note that this reduces to the noise correlation function for the equilibrium bulk noise when 
k± » 1/d, i.e. when the membrane fluctuations are not significantly affected by the wall 
(see Section 111). 

We note here that the equation of motion of the height field can be derived in the limit 
1/L C << k± « 1/d by directly imposing impermeability as in [29]. Since the membrane is 
impermeable, the height field is conserved and obeys a continuity equation d t h = —V±.jh- 
The height "current" can be calculated from the lubrication approximation as jh = Vj_P, 
where P is the pressure due to the elastic forces and the noise. Using the appropriate 
expressions for these terms, we recover the equation for the height (Eq. (4.4)) in the "close 
to the wall" limit. 

The eigen-modes follow from Eq. (4.6) 

u u u 2 = ^ (Dk 2 ± + °<3kl ± K gki) + X - [\j{cyQk\ + Dk\ + K<3k\f - AQD(a + K e k 2 ± )kij 

In the absence of curvature-concentration coupling the modes corresponding to the hydro- 
dynamic relaxation of a two-component membrane near a wall are recovered : Wi,U2 — 
Q{akl + Kkj_), Dkl- For membranes with curvature-concentration coupling the modes are 
modified to 

ui w Q{okl + K e k\),u 2 w Dkl + G^—zrkl. (4.10) 

Xo 

It is easy to see that these reduce to the expressions derived in the bulk case, in the limit in 
which membrane fluctuations are unaffected by the wall. 

V. NON-EQUILIBRIUM MEMBRANES: MODE STRUCTURE IN BULK 

We now discuss the dynamic correlations and mode structure of an impermeable two- 
component membrane with active force-centres, when placed far away from confining walls. 
For inactive force centres, the static properties of such a system are given by the free energy 
of Eq. (2.2). However, in the presence of non-equilibrium driving forces, thermodynamic 
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restoring forces must be supplemented by additional terms. These are discussed briefly 
here. 

First, we must account for non-equilibrium forces. In the context of the membrane-pump 
system, these arise due to the action of the pumps. However, since pumping refers to forces 
internal to the system, the force density exerted on the whole system (membrane + pump 
+ solvent) must integrate to zero on length scales comparable to the pump size. A simple 
representation of this requirement idealizes the pump as two force-centres of opposite sign 
but equal magnitude separated by a distance + and embedded in the bilayer. Such 
a "force dipole" representation was introduced in Ref. [28] and will represent an important 
part of our discussion here. 

Second, we must account for the non-equilibrium nature of the noise present in the fluid 
due to the switching "on" or "off" of the pumps. For simplicity, we assume that this noise 
is delta correlated. Such an assumption corresponds to a coarse-graining in time across 
the typical pump "dead" time which separates two pumping events. The non-equilibrium 
active noise can now be absorbed into the correlation function of the equilibrium noise 
(Eq.(3.9)). However, note that the "temperature" which appears in this correlator is not 
the thermodynamic temperature. To smoothly interpolate to the equilibrium case, we must 
redefine this temperature as an effective temperature which becomes the true thermodynamic 
temperature in the absence of activity. 

Several of the results in this section overlap with those of Manneville et al. (Ref. [28]) in 
which the effects of permeation and curvature induced activity on the dynamics of an active 
membrane in bulk are considered. In our model we assume impermeability from the outset 
and neglect the effects of curvature induced activity. Such a "minimal" model makes the 
calculation of the correlation functions of a wall-bounded active membrane more tractable 
and transparent. 

We begin by writing the equations of motion for the various fields. The fluid 
is incompressible and hence obeys Eq. (3.10). Since no external forces act on the 
(fluid+membrane+proteins) system, momentum is conserved [42,44] and Eq.(3.7) holds. 
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The momentum current tensor tvh is then 



■Kij = P5ij + pViVj - //(V,/' ; + VjVi) + S i:j + A, 



(5.1) 



where v(r, t) is the three dimensional fluid velocity, P(r, t) is the three dimensional pressure 
field. Sij is the contribution to the stress tensor due to membrane conformations and is 
the contribution from the active pumps. Eqns. (3.8) and (3.10) can be solved as before for 
the transverse (£), longitudinal (/) and z components of the fluid velocity, in terms of the 
t,l and z components of the force F(r, t) acting on the fluid due to membrane and pump 
stresses. The transverse components will be ignored as before; fh(r, t) represents the thermal 
noise present in the bulk of the fluid. 

We can now calculate v z (k±,0, t) by solving the equations of motion with appropriate 
boundary conditions. The equation of motion for the height field in the absence of perme- 
ation is d t h(r±,t) = v z (r±,0,t), where v z (r±,0,t) is the velocity of the surrounding fluid at 
the mean position of the membrane. The membrane does not feel the forces from the pumps 
directly because of the absence of permeation but only via the forces the pumps exert on 
the fluid [28,29]. The contribution of the bulk equilibrium thermal noise present in the fluid 
to the z— component of the fluid velocity at z = is f z . The statistics of f z have been 
discussed earlier (see Section III). The density difference field obeys the diffusion equation 
Eq. (3.2). As in the equilibrium case, we work in the linearly stable regime. 

The membrane fluctuates about z — 0. The force exerted by the membrane and pumps 
on the fluid is of the form 



The force f pU m P -fiuid is modelled as a dipolar force density with force-centres located asym- 
metrically about the z = plane. 



F(k ± ,/) = ^)[ 



5_F_ 

Jh 



Z + Fi kj_] + fpump-fluidT 1 - 



(5.2) 



fpump-fiuid = fip(r±,t)[5(z' - w 



1)-5(z' + w% 



(5.3) 



where / sets the scale for the active force. 
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We again solve for the unknown coefficient F™ using membrane incompressibility. This 
imposes vi(k±,0,t) = in Eq. (3.11). The hydro dynamic velocity obtained at z = using 
Eq. (3.11) enters in the equation of motion for the height field: 



d t h(k ± ,t) 



1 r) P 



(5.4) 



where 



Q(k ± ,w\w l ) = -e- k±w \l + k ± w l ) + e- k±wt (l + k ± w^ 

= (itT -W) JH (UT -W ) + .... 



(5.5) 
(5.6) 



Note that Q(k±, ,w^) vanishes when = u^, implying that an asymmetry in the position 
of the force-centres above and below the membrane mid-plane is needed for a non-zero active 
(non-equilibrium) force. The equations of motion are 



1 h(k ± ,u)^ 



4r/k ± 

\KHk 4 t 



4r)k ± 



\ 


' h(k ± ,u)^ 




' /,(kj 








+ 




J 











(5.7) 



The correlation function for the height fluctuations of an active membrane (about the flat, 
stable, steady state phase) is then 

i , pn* x 

« ^ 4r,k ± Xo Lj 



(h*(k ± ,t)h(k ± ,t)) 



k R T ir^ + r h 1 -T- 1 + -^ 



(5.8) 



W ± (V 1 + 0(r ft - ie + r-i) » 

where the quantities r^ 1 , r^ 1 have been defined earlier. r^ ie has the same form as r^ 1 with 
k replaced by K e . The relaxation time 



1 kH 



T fQ(k ± ,w\w l )k 2 



(5.9) 



We discuss the positivity of the right hand-side of Eq. (5.8) after discussing the mode 
structure of the system. From the equal time correlation function we can calculate the 
quantity (h 2 (L)). To lowest order we find 

kHVw' 



(h\L)) = k B TL' 



Xo' 1 



(5.10) 
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for tensionless membranes and 



(h\L)) = ^\n(L/a), (5.11) 
a 



for tense membranes. V = ^ w \~ wl \ where w is the size of the pump, represents the force 
dipole energy. It is estimated to be a few ksT [28]. These results imply that on large 
length scales, pumping activity does not change the roughness of equilibrium, tensionless 
membranes but shifts the temperature at lowest order to 

kHVw 



1 + 



(5.12) 



For tense membranes both roughness and effective temperature are not affected by non- 
equilibrium effects at large length scales since the divergence of the height correlation func- 
tion with system size remains logarithmic. 

In the stable regime Eq. (5.10) is positive, as discussed later in this section. The correla- 
tion function in Eq. (5.10) is calculated from the lowest order diffusive terms. Manneville et 
al. [28] calculate an identical correlation function but neglect these diffusive terms in favor 
of higher order terms on the grounds that the diffusion constant is small in the physical 
situation. For this reason the effective temperature we calculate in Eq. (5.12) is not the 
effective temperature that Manneville et al. calculate and compare with experiment. 

The height correlation functions we calculate in our linearized approach (Eqns. (5.10) 
and (5.11)) show that the roughness of an active fluid membrane is unchanged in the bulk. 
The calculation of Ref. [29] obtains an enhanced roughening. The reason is the different 
nature of the force densities used in both the calculations. Prost et al. consider a model in 
which the membrane feels the activity of the pump only as a consequence of permeation. In 
our calculation, Eq. (5.4) written as a Langevin equation is 

d t h(k±,t) + r^h(k ± ,t) = - ^"^V + U (5-13) 

We assume that the pump is always in the active state over the time scales of our interest. 
The forces on the right-hand side of the above equation arise from three sources: an elastic 
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part (due to the curvature concentration-coupling term in the free energy), a asymmetric 
dipole part (from the activity of the pumps) and the thermal noise from the fluid. Since 
the membrane is impermeable it does not feel the active forces or thermal noise directly 
but only via the fluid. We emphasize that in the present model the membrane senses active 
forces and thermal noise even when \ p = 0, in contrast to the model of Prost et. al. 
The modes are 

(jj\,uj2 = 3— (ck± + ArjDk\ + Kk± 3 



877 

± (V (° k ± + 4f]Dk\ + Kk ± 3 ) 2 - 16r](Dn e kl + Dak\ + fUnMlk])^ . 

In the absence of curvature-concentration coupling, non-equilibrium effects vanish to lowest 
order and the modes are those for an equilibrium membrane with hydrodynamic interactions 
ui,u>2 = (crk± + K,k\) / 477, Dk\. When a — 0, curvature-concentration coupling modifies this 
mode structure to 

K e , o fKH(w^ 2 — W^ 2 ) , o 

„ kH(2kH - f(w l2 -w^))k\ , ri<( , 

uj 2 « Dkl + ^ -P — 5.14 

877x0 1 

We see that to lowest order, non-equilibrium effects are apparent only at non-zero curvature- 
concentration coupling. We also observe that the non-equilibrium terms vanish when = 

For stability we require that all the eigenmodes of the system must decay. This means 

that uji and uo 2 are positive. To the order above, this imposes ^ — 1 < < ^. If HV is 

K x K 

such that this condition is satisfied then the system is stable independent of the sign of this 
parameter. Typical experimental values such as k = 10 fc^T, kH = w /cbT,Xo 1 = w 2 ksT, 
f = 10"~ 12 iV, satisfy the above conditions [39]. 
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VI. NON-EQUILIBRIUM MEMBRANES : MODE STRUCTURE IN THE 

PRESENCE OF A WALL 



We now study the effect of a wall located at z = on the fluctuations of an active 
membrane located at z — d. The forces exerted on the fluid by the membrane, pumps and 
the wall are now given by 

SF 

F{k ± ,z') = 5{z' - d)[-—z + F/"kJ + fif>(T ± ,t)[6(z' -d-w v )- 5{z' -d + w l )]z 

+ 5(z / )[F?z + F l w k ± }. (6.1) 

We need to determine the unknown coefficients F( n ) Ff and F™ . The following three bound- 
ary conditions are used, following Eqns. (3.11), 

v z (k ± , 0, t) = v t (k ± , 0, t) = v t (k ± , d, t) = 0. (6.2) 

These boundary conditions impose the following constraints: (i) v z (k±,z = 0) = 0, i.e. the 
^-component of the velocity vanishes at the wall (ii) vi(k±, z = 0) = 0, i.e. the /-component 
of the velocity vanishes at the wall, and (iii) vi(k±,z = d) = 0, i.e. the membrane located 
at z — d is incompressible. We neglect the tangential component as before. 

The values of the three unknown coefficients are obtained on solving the three simulta- 
neous equations (Eq. (3.11)) together with the three boundary conditions. These values 
can then be used to obtain the velocity of the fluid at the location of the membrane. The 
velocity thus obtained has a part arising from equilibrium elastic forces and a part arising 
from non-equilibrium active forces. These can be separated out, yielding the equation of 
motion for the height field as 

d t h(k ± ,t) = v z (k ± ,d,t) 

= -g( 6 -^) + g aw (N) + f z (k ± ,t), (6.3) 

where Q is the usual hydrodynamic kernel obtained for the equilibrium fluctuations in the 
presence of a wall and the new non- equilibrium term 
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g aw (k ± ,d,w\w l ) 



j^(3dw l - 2w l ) + 0(ki) + ... if k ± < 

if A; ± > l/d. 



(6.4) 



is the hydrodynamic kernel for non-equilibrium fluctuations in the presence of the wall 
[29] (the expanded form of this kernel is given in [45]; see also Ref. [46]). In the expressions 
for Q and Q aw , the lower cut-off on k± is the inverse of the collision length L C) which we 
will calculate towards the end of this section. The membrane-wall distances we use in our 
calculations are such that d » w\w^ and hence w^ 3 can be neglected in comparison to 
dw^- 2 . In the long wavelength limit, Q aw depends on and not on (w^ —w 1 *) as in the bulk. 
When the fluctuations are not affected by the wall, Q aw vanishes when equals u4 since 
in bulk the directions "up" and "down" are equivalent. The non-equilibrium contribution 
to the force density therefore vanishes. Thus, in this large membrane- wall distance limit, we 
recover the previous situation of active membrane fluctuations where the asymmetry was 
crucial for non-zero active force density. The wall breaks the symmetry in the "up" and 
"down" directions and hence the hydrodynamic kernel picks out a particular direction and 
does not depend on (w^ — w^) when membrane fluctuations feel the wall. The statistics for 
the noise f z in the presence of the wall follows from Eqns. (4.8). 
The equal time correlation function can be calculated from 



h(k ± ,u] 
^(k_L,w) 



( i 



+ 



ik i .£ 



(6.5) 



4> J 



The correlation function is 



(h*(k ± )h(k ± )) = k B TQ 



iw f 2 g aw2 ^\ 



(6.6) 



where 



Q aw —Jk\. 

Xo 



(6.7) 



r h ' has the same form as r h with k replaced by K e . Inserting physically reasonable 
values for the parameters which enter these definitions (see below), we verify that the right 
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hand-side of the correlation function is always positive. Thus, this correlation function is 
calculated in a regime in which the active membrane is in a stable steady-state phase. 
At large length scales and in the presence of a confining wall (1/L C <C kj_ <C l/d), 



cf D 



ln(L/a) 



(6.8) 



for tensionless membranes and 

e' + ifeyDxo- 1 ] 



(h\L)) = k B T 



k B T 

ln(L/a) ~ ln(L/a), 



2 



(6.9) 



e'a + iefKH/xo- 1 ), 

for tense membranes. We have defined the quantities e = dw lZ /(4i]) and e' = d 3 /(l2i]) for 
notational convenience. Thus, the presence of the wall makes tensionless active membranes 
less rough than they are in bulk since the height correlation function scales now logarithmi- 
cally with L instead of quadratically. Note that this effect is opposite to that predicted via 
the scaling arguments of Prost et al. in Ref. [29]. 
The effective surface tension is then 



kH 



e'xo' 1 fe 
ef D 



(6.10) 



This is an explicitly non-equilibrium effect, since E = when / = 0. The induced surface 
tension also crucially depends on the curvature-concentration coupling H since it vanishes 
on setting this coupling to zero. Therefore, both non-equilibrium forces and curvature- 
concentration coupling are essential for the stiffening of a tensionless, active membrane near 
a wall. 

To calculate the collision length L c we set h 2 (L c ) pa d 2 in Eqns. (6.8) and (6.9). Thus, 



L c (d) 



a exp {Y,d 2 /k B T) if o = 0, 
a exp (a e d 2 /k B T) if a ^ 0. 

The eigen-modes have frequencies 



(6.11) 



2 

± i ( \l [Dk] -f Qak\ ■■■■ g,-d, [Y : - l{fg""'HH\k[ GDh,^ GDal A _) 
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In the absence of curvature-concentration coupling, non-equilibrium forces have no effect on 
the modes and the modes reduce to those of an equilibrium membrane in the presence of 
a wall, within the linearized theory. For tensionless membranes with non-zero curvature- 
concentration coupling, non-equilibrium effects become apparent in the mode structure. The 
modes to lowest order are 

Xo 

^a;,**!^^ (6 . 12) 

Xo 

When membrane fluctuations are influenced by the wall, Q = d 3 k\/12r] and Q aw = 
dw^k^/An, to lowest order. Far away from the wall the modes reduce to those of a non- 
equilibrium membrane in the bulk. 

It may seem that for initially tensionless membranes near a wall, the induced surface ten- 
sion E can have either sign depending on the sign of the curvature concentration coupling 
H. We will show that for stability, it is imperative that H is positive for active membranes 
close to walls; a negative H leads to a negative u>i, and thus an instability. The require- 
ment of positivity of H in the stable regime ensures that our approximations of considering 
fluctuations at linear order about an initial flat steady-state phase are self-consistent. 

We now present a physical argument for why pumping activity should lead to the stiff- 
ening of active, tensionless membranes near repulsive walls. The calculations of Section IV 
indicate that the eigen-values governing the relaxation of a tense, equilibrium membrane 
near a wall are ad 3 k^_/(12i]) and Dk\. For an active, tensionless membrane near a wall 
these eigenvalues are KH f dw k ± anc [ _ rj^g coe fg c j eri t Q f the slower mode for the active 
case (0(k\)) can be written as ^^■d 3 /(12i]) and compared to the expression for the slow- 
est relaxation mode for tense, equilibrium membranes. This comparison suggests that an 
apparent surface tension ( K f H j{ ) is induced due to purely non-equilibrium effects, even for a 
membrane with vanishing bare surface tension. 

In the linearized theory, the induced surface tension obtained using the above argument 
and the surface tension E obtained from the coefficient of the logarithmic term in the cal- 
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culation of the height roughness coincide. It is thus apparent that the presence of the wall 
smoothens a tensionless, active membrane by making it relax faster than its tensionless 
equilibrium counterpart (whose slowest mode is 0(kj_)). 

Why does an active membrane relax faster near the wall? Height field fluctuations are 
predominantly governed by the tp field (the second term in Eq. 6.3) because this term is 
0(k"j_) in contrast to the other terms which are of higher order in wavenumber. Thus, this 
equation to lowest order in wavenumber is d t h = G aw f^>- Since the ip field relaxation is fast 
compared to the height field, we may assume that, over the larger time-scales of interest to 
us in studying the dynamics of the height field, ip relaxes instantaneously to ip s = — ^^rk^h. 
The effective height field equation now becomes 

d t h = -^ r k A ± h. (6.13) 

Xo 

Note that the coefficient of the k\ term in the above equation can be rewritten to obtain 
the induced surface tension E, as explained in the previous paragraph. 

The above argument invokes and crucially requires the presence of the wall, activity and 
curvature-concentration coupling. In physical terms, height field fluctuations are accompa- 
nied by fast relaxation of the protein density field ip. This relaxation is to a steady state 
value determined by the pump diffusion coefficient and curvature-concentration coupling. 
The pumping activity is then inhomogeneously distributed along the membrane surface. 
To see the physical consequences of this inhomogeneous distribution, consider a fluctuation 
creating a region of positive curvature. If H is positive, the number of "up" pumps in the 
region of positive curvature tends to decrease with time and the decrease in the pumping 
activity causes a smoothening of the membrane. However, if H is negative, the region of pos- 
itive curvature attracts more "up" pumps which further bootstraps the height fluctuation, 
leading to an instability. 

The explicit expressions obtained above also permit us to estimate the collision length 
for tensionless active membranes in the presence of repulsive walls. Using estimates given 
in [28] for various physical quantities (re = 10~ 20 Joules, H = 5 nm, f = 1CT 12 N,d — 
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1 micron, xo 1 = 1CT 21 Joule m 2 ), we find that the induced surface tension £ ~ 1CT 7 — 
1CT 8 N/m. This translates to collision lengths L c of the order of 10 5 — 10 6 times molecular 
scales. Assuming molecular length scales of order a nanometer, this calculation would then 
predict collision lengths of the order of millimeters, for membranes with physical parameters 
in the above range. 

VII. CONCLUSIONS 

This paper has presented a calculation of the linearized hydrodynamics of two-component 
membranes including the effects of curvature-concentration coupling, specifically in the con- 
text of a simple model for a biological membrane with protein and lipid constituents. We 
have studied in some detail the case in which such membranes are placed close to as well as 
far away from a confining repulsive wall. We have considered the case in which such proteins 
are "active", in which case they exert non-equilibrium forces on the solvent, as well as the 
case of "inactive" pumps, corresponding to thermal equilibrium. 

The calculations we report are for a model which shares many similarities with that of 
Ref. [28]. Manneville et al. consider the more general case of a membrane with a finite 
permeability whereas we work with an impermeable membrane throughout. In addition, 
we neglect the possible effects of curvature on activity. However, the relatively simpler 
structure of our model enables us to establish a connection with equilibrium results in 
a more transparent fashion. We work in a framework in which boundary conditions are 
explicitly implemented. This is especially of use when we consider the linearized theory of 
fluctuations of an active membrane near a wall. The results we report for the linearized 
hydrodynamics of active and inactive two-component membranes near walls are new, as are 
our results for the equilibrium case far away from confining walls. 

We find, far away from the walls and consistent with the results of Manneville et al. [28], 
that the presence of active pumps on the membrane does not change its roughness. In the 
case of tensionless membranes, fluctuations can be described in terms of an effective tem- 
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perature which is shifted from the thermodynamic temperature by non-equilibrium effects 
due to pumping. These terms depend on the quantity (w^ — w 1 *), a measure of the degree of 
asymmetry in the positions of the force-centres. We reiterate the observation of Manneville 
et al. in Ref. [28] that this asymmetry is essential for non-equilibrium effects to be visible 
in the membrane dynamics. 

We find that the activity of the pumps within the stable regime appears to smoothen out 
height fluctuations in the presence of the wall since height correlations scale only logarith- 
mically with L for both tense and tensionless membranes. This smoothening out within the 
linear theory presumably competes with the possible enhancement of fluctuations, if any, 
which might arise in a theory which incorporates the effects of non-linearities. We have 
provided a physical interpretation of this effect in this paper and obtained estimates for the 
appropriate collision length. We find that this collision length can be significantly larger than 
typical length scales for biological vesicles, indicating the relevance of the physical effects 
we study here to a complete description of the hydrodynamics of active pump-membrane 
systems. It is important to note that this smoothening relies crucially on the existence of 
curvature-concentration coupling terms, neglected in previous work on the near-wall case. 
These results suggest that earlier ones based on scaling arguments which predicted an en- 
hancement in the roughness may need to be reevaluated in the light of this work. 

We have also calculated the "crossover" function Q aw and linearized mode structure in the 
presence of the wall. While the positional asymmetry in the placement of force-centres was 
crucial for active membrane fluctuations in bulk, we point out that non-equilibrium effects 
are apparent even for symmetrically placed pumps in the presence of a wall. This unusual, 
although not totally unexpected, result does not appear to have been noticed before. 

We also find that the requirement of stability for active membranes near a wall constrains 
H, the coefficient coupling concentration with curvature, to be positive. This can have 
interesting consequences. Consider two active membranes, one with positive H and the 
other with negative H, fluctuating in the bulk. Provided HV for each lies in the window 
of stability in the bulk, these membranes are stable. Now bring both close to a wall. The 

31 



membrane with positive H will remain stable. The fluctuations of the membrane with 
negative H however, will turn unstable; its final fate will be decided by the non-linear terms 
we have omitted from our discussion. 

Because our calculation is an explicitly linear one, we cannot address issues such as the 
relevance of such non-linear terms. The arguments of Prost, Manneville and Bruinsma which 
predict an anomalous roughening of active membranes near walls rely on identifying the role 
of such non-linear terms and estimations of their importance via scaling arguments. It 
would be interesting to construct similar arguments for the model studied here, particularly 
since our results in the linear case indicate that fluctuations of an active membrane near 
a repulsive wall can be smoothened in some regimes, rather than roughened relative to the 
equilibrium case. Further work along these lines is currently in progress. 
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VIII. APPENDIX 



A. A note on the impermeability of the membrane 

In this sub-section we present an argument, following Manneville et. al, to justify the 
assumption of membrane impermeability. Had we incorporated permeation due to activity 
in our calculations, active-permeative and active- hydrodynamic terms would appear as (A p — 



' Sri W i\ ls the permeation coefficient) [28]. The cross-over length scale from the 



hydrodynamic to the permeative regime occurs when k± = w( -^t_ p w I) ; *- e - l± = nw( -^~ wl ^ (w 
is the size of the pump). With A p = 10" 12 m 3 /N.s,rj = 10~ 3 kg/m.s and w = 5.1CT 9 m [28], 
this crossover length scale is estimated to be of the order of a centimeter. The length scales 
of concern to us are in the micron and sub- micron range [28]. Thus, at these length scales, 
effects due to permeation can be ignored. These arguments can be extended to estimating 
the role of passive permeative vs. hydrodynamic modes of relaxation; we conclude that 
at the length scales of interest to us, the assumption of membrane impermeability remains 
valid. 

In the presence of a wall, active permeative terms would dominate if \ p > Q aw . This 



/ d I 2 

means that the length scale above which permeation dominates is J n 2 ™ x . This length is 
approximately 40 microns when the membrane-wall distance d = 1 micron. Hence, in the 
regimes of interest to us, permeation can be ignored and relaxation is purely hydrodynamic. 
Hence, we assume \ = to derive our results. 

B. A modified model for active membranes 

We present here results for a model in which proteins are located asymmetrically across 
the bilayer with their centers of mass displaced above or below the bilayer midpoint (see 
Fig. 3). The head and tail of the pumps are then located at distances W[ and w s respectively 
from the midpoint of the bilayer. Note that in this model, the "up" and "down" pumps are 
related to each other by reflection along the z = line unlike the model of Fig. 2. This 
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model can be used to describe situations in which proteins are attached to planar bilayer 
membranes which are symmetric on both sides, unlike vesicles in which the inner leaf and 
outer leaf need not be symmetrically preferred by the protein. 

When pumps are active on such a membrane in bulk, the force density due to the pumps 
acting on the fluid has the form 

fpump-fiuid = f[6(z - wi)n} - 5(z + w s )n T - S(z + wi)n l + S(z - w s )n l ], (8.1) 

where w) and are the local number of "up" and "down" pumps respectively. On re- 
peating the calculations done in Section V, we find that the z-component of the hydro- 
dynamic velocity has the same form as shown in Eq. (5.7) with Q(k±,w\w^) replaced by 
tt'(k± : w\, w s ) = —e k±Wa (l + k±w s ) + e k±Wl (l + k±wi), which vanishes when wi = w s (i.e. for 
symmetrically placed proteins). Correlation functions and relaxation times can be read off 
from the corresponding expressions in Section V by identifying and w 1 * with w s and wi 
respectively. For example, 

/f T -i | r -i r -i | Z 2 " 2 A 



(h*(k ± ,t)h(k ± ,t)) = ^ 



(8.2) 

y \' i/j ' ' ft J \ ' ft ' ' a J J 

where the quantities r^ 1 , r^ 1 have been defined earlier, r^ 1 has the same form as r^ 1 with 
n replaced by K e . The relaxation time 

1 K~H 

r- 1 = /flfo,, (8.3) 

We note here that this identification is purely notational since (w^w 1 *) and (w s ,wi) physi- 
cally represent different lengths in the two models. 

When the calculations are repeated for such a model in the presence of a repulsive 
wall, we find that the hydrodynamic velocity acquires an additional contribution from the <fi 
(total protein density) field. In this linearized calculation we assume that the compressibility 
associated with the field is very large and ignore its dynamics. However, we do not expect 
this contribution to alter the long wavelength fluctuations of the active membrane near the 
wall significantly. 
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Figures 



z=0 



Fig. 1. A schematic diagram of the two kinds of proteins, represented as "up" and "down" 
triangles. These affect the local curvature of the membrane by bending it towards or away from 
the local normal. The bilayer mid-plane lies at z = 0. 



Fig. 2. The asymmetric dipole model used to describe the activity of the proteins. The 
force-centres are places at distances and from the bilayer midpoint. The centers of mass of 
the pumps are displaced above the bilayer midpoint. The superposed "up" and "down" triangles 
indicate the underlying shape asymmetry of the two different kinds of pumps. 
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Fig. 3. The asymmetric dipole model of the pumps described in the Appendix. The centers of 
mass of the pumps are displaced above or below the bilayer midpoint. The distance of the furthest 
force-centre from z = is wi and the distance of the nearest force-centre from z = is w s . 
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